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SUMMARY 


Sufficient  conditions  are  obtained  which  provide  that  a lifetime 
density  has  a bathtub  shaped  failure  rate.  Analogous  conditions  handle 
increasing,  decreasing,  and  upside-down  bathtub  shaped  failure  rates. 
Application  of  these  results  to  exponential  families  of  densities  is 
particularly  straightforward  and  effective.  Examples  are  furnished 
which  introduce  new  bathtub  models  and  illustrate  the  use  of  the  general 
results  for  existing  models.  Examples  involving  mixtures  are  considered. 
Maximum  likelihood  estimation  for  one  of  the  bathtub  models  is  described. 

1 . INTRODUCTION 

The  probability  distribution  of  the  time-to-f ailure  of  an  item  can 
be  characterized  by  the  failure  rate,  r(t)  = f(t)/R(t),  where  f(t)  denotes 
the  density  function  and  R(t)  the  reliability,  or  probability  of  failure 
after  time  t.  The  failure  rate  has  a probabilistic  interpretation: 
r(t)dt  represents  the  probability  that  an  item  of  age  t will  fail  in 
the  interval  (t,  t + dt). 

Many  parametric  lifetime  models,  such  as  the  gamma,  Weibull,  and 
truncated  normal  distributions,  have  monotone  failure  rate.  If  r(t) 
increases  monotonically  over  time,  the  distribution  is  said  to  have 
increasing  failure  rate  (IFR).  If  r(t)  decreases  monotonically,  we 
have  decreasing  failure  rate  (DFR).  The  IFR  property  is  characteristic 
of  devices  which  consistently  deteriorate  with  age,  whereas  the  DFR 
property  is  characteristic  of  devices  which  consistently  improve  with 
age.  Many  physical  phenomena  exhibit  failure  rates  which  are  non- 
monotonic. A common  description,  which  is  appropriate  for  modeling 
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human  lifetimes,  shows  three  phases:  an  initial  phase  where  the  failure 
rate  decreases,  followed  by  a middle  phase  where  the  failure  rate  is 
essentially  constant,  concluded  by  a final  phase  where  the  failure  rate 
increases.  For  humans,  the  first  (infant  mortality)  phase  shows  deaths 
due  to  hereditary  defects,  whose  impact  dimishes  with  time.  The  middle 
(chance  failure)  phase  shows  deaths  due  typically  to  sudden  jolts,  such 
as  accidents.  The  final  (wear-out)  phase  shows  death  resulting  from 
the  natural  accumulation  of  negative  effects.  Such  failure  rates  are 
usually  termed  bathtub  (BT)  shaped.  The  logical  counterpart  to  BT 
failure  rate  is  the  three  phase  situation  where  the  failure  rate  initially 
increases,  then  becomes  essentially  constant,  and  ultimately  decreases. 
This  failure  rate  function,  which  we  will  term  upside-down  bathtub  (UBT) 
shaped,  can  be  found  in  accelerated  life  testing,  where  the  items  tested 
are  subjected  to  abnormally  high  stress  levels.  The  lognormal  and  inverse 
Gaussian  lifetime  models,  as  we  shall  show,  have  UBT  shaped  failure  rates. 
For  the  sake  of  brevity  in  this  paper,  we  shall  frequently  say  that  the 
failure  time  model,  distribution,  or  f is  BT,  IFR,  UBT,  or  DFR,  when 
we  mean  more  precisely  that  the  associated  failure  rate  function  is  BT 
shaped,  increasing,  etc. 

Although  bathtub  shaped  lifetime  models  are  of  great  practical 
value,  few  have  been  suggested  in  the  literature.  See  Lieberman  (1969). 
One  reason  for  this  scarcity  undoubtedly  has  been  the  difficulty  in 
ascertaining  whether  a given  r(t)  is  bathtub  shaped.  So  motivated,  we 
obtain  in  Section  2 a general  result  which  furnishes  sufficient  con- 
ditions that  a distribution  has  a BT  shaped  failure  rate  function.  Dual 
results  handle  the  UBT,  IFR,  and  DFR  situations.  For  lifetime  distribu- 
tions of  the  exponential  family  type,  or  mixtures  of  the  same,  these 
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conditions  have  a simple,  easy  to  test  form.  This  ease  of  application 
is  demonstrated  in  Section  3,  where  new  BT  lifetime  models  are  proposed. 

In  Section  4,  statistical  inference  is  considered  for  the  new  models. 

The  applicability  of  the  general  results  to  mixtures  is  useful. 

In  certain  situations,  it  is  reasonable  to  assume  that  a proportion  p 
of  items  in  a population  come  from  one  distribution  (say  IFR,  characterized 
by  wear-out  failures),  and  the  remaining  proportion,  q = 1 - p,  of  items 
come  from  another  distribution  (say  DFR,  characterized  by  catastrophic 
failures).  An  example  with  electron  tubes  is  considered  by  Kao  (1959). 

We  see  in  Section  3 that  for  certain  gamma  mixtures  the  combined  popula- 
tion has  a BT  shaped  failure  rate  function. 

We  shall  assume  throughout  this  paper  that  the  failure  time  dis- 
tribution has  support  (O.M),  i.e.,  (t:f(t)  > 0}  = (0,M).  The  constant 

M may  be,  and  indeed  in  most  applications  will  be,  «.  Because  of  further 

assumptions  placed  on  the  density  function  in  Section  2,  the  failure  rate 
function  r(t)  will  be  continuous  and  twice  differentiable  for  all 
t € (0,M).  With  exception  of  the  exponential  distribution,  where 
r' (t)  = 0 for  all  t > 0,  resulting  failure  rate  functions  will  exhibit 
properties  of  strict  monotonicity.  For  brevity,  we  denote  by  (I)  the 
strict  IFR  situation,  "r1 (t)  > 0 for  all  t€  (0,M),"  and  denote  by  (D) 
the  strict  DFR  situation,  "r' (t)  < 0 for  all  t€  (0,M)."  Similarly, 
we  denoto  by  (B)  the  special  BT  situation,  "for  some  t* € (0,M),  r' (t)  < 0 
for  all  t € (0, t*) , r' (t*)  * 0,  and  r' (t)  > 0 for  all  t € (t*,M)." 

Finally,  we  denote  by  (U)  the  special  UBT  situation,  "for  some  t*€  (0,M), 

r' (t)  > 0 for  all  t € (0,t*),  r' (t*)  = 0,  and  r' (t)  < 0 for  all 

t € (t*,M) ." 


4 


2.  GENERAL  RESULTS 

In  this  section  we  shall  obtain  general  results  which  supply 

sufficient  conditions  to  characterize  a given  failure  time  distribution 

as  being  either  BT,  IFR,  UBT,  or  DFR.  We  assume  throughout  that  the 

failure  time  distribution  is  absolutely  continuous  with  support  (0,M), 

where  the  constant  M may  be  «.  We  assume  the  density  f(t)  is  con- 

•k 

tinous  and  twice  differentiable  on  (0,M).  We  define  g(t)  as  the 
reciprocal  of  the  failure  rate, 

g(t)  = l/r(t)  = R(t)/f(t).  (2.1) 

It  follows  that  g(t)  is  positive  valued,  continuous,  and  twice  differen- 
tiable on  (0,M).  In  fact,  we  have 

g’(t)  - g(t)7](t)  - 1,  (2.2) 

where  T](t)  is  defined  by 

71(t)  = - f'(t)/f(t).  (2.3) 

THEOREM.  Suppose  g' (t)  can  be  expressed  as 

,M 

g'(t)  - s (y)  [T)(t ) - Tl(y)]dy,  (2.4) 

“t 

for  some  function  st(y)  which  is  positive  valued  for  all  t and  y in 
(0,M). 

(a)  If  T|'  (t)  > 0 for  all  t € (0,M),  then  (I). 

(b)  If  Tl'(t)  < 0 for  all  t€  (0,M),  then  (D). 

♦All  differentiation  mentioned  and  performed  in  this  section  will  refer  to 
partial  differentiation  with  respect  to  the  time  argument  t;  i.e.,  all  para- 
meter values  are  held  fixed.  Such  derivatives  will  be  denoted  by  ' and  ". 
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(c) 


and 


Suppose  there  exists  tg€(0,M)  such  that 
tl'(t)  < 0 for  all  t€  <0,t0),  TJ' (tQ)  = 0, 
Tl*  (t)  > 0 for  all  t€  (tg,M). 


(2.5) 


(i)  If  there  exists  yg€(0,M)  such  that  g'(yg)  = 0,  then  (B). 

(ii)  If  there  does  not  exist  yQ€  (0,M)  such  that  g'(yg)=0,  then  (I). 


(d) 


and 


Suppose 

Tl'(t)  > 

T| ' (t)  < 

(i)  If 
(ii)  If 


there  exists  tg€(0,M)  such  that 
0 for  all  t€  (0,t0),  T)'  (tQ)  = 0, 

(2.6) 

0 for  all  t € (tg,M) . 

there  exists  y^€  (0,M)  such  that  g'(yg)=0,  then  (U). 
there  does  not  exist  y^S  (0,M)  such  that  g'(yg)=0,  then  (D). 


PROOF.  (a)  The  assumption  that  T|'  (t)  > 0 for  all  tc  (0,M)  implies,  from 
(2.4),  that  g'(t)  < 0 for  all  t€(0,M),  which,  from  (2.1),  implies  (I). 

(b)  Here  it  follows  that  g' (t)  > 0 for  all  t€(0,M),  which 
implies  (D). 

(c,i)  Claim  g"(yg)  < 0.  Since  g'(yg)  = 0,  it  follows  from  dif- 
ferentiation of  (2.2)  that  g"(y0)  = g(y0)‘n' (yQ)  • Therefore  g"(y0)  < 0 ® 
T|*  (yg)  < 0 ® y^  < tQ.  Suppose  y g > tQ.  From  (2.4)  it  is  apparent  that 
g'(t)  < 0 for  all  t€[tg,M).  Therefore  g'(yg)  < 0,  which  is  a contra- 
diction. Hence  y g < Cg  and  g"(yg)  < 0.  I»-  is  clear  that  there  is  only 

one  root  in  (0,M)  to  g'(y)  = 0,  namely  y=yg,  and  g attains  a maximum 
at  this  point.  This  implies  (B),  with  t*  = yQ. 
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(c,ii)  Here  either  g'(t)>  0 for  all  t€(0,M)or  g'(t)  < 0 for  all 
t€(0,M).  From  (2.4)  we  have  that  g'(t)  < 0 for  all  t€[tQ,M).  Therefore 
g'(t)  < 0 for  all  t€(0,M),  and  (I)  holds. 

(d)  The  proof  is  analogous  to  that  of  (c)  and  will  be  omitted. 


The  Theorem  is  readily  applicable  to  exponential  families  of 
densities.  Suppose  f can  be  expressed  as 


f(t)  = C (9 ) expj^=10.Ui(t)J,  0 < t < M, 


(2.7) 


where  for  each  i,  IL(t)  is  twice  differentiable  on  (0,M).  Then  from 


(2.3), 


^(t)  = ■^10iui(t).  and  Tl'(t)- -2i=1eiui(t). 


(2.8) 


The  condition  (2.4)  of  the  theorem  always  holds  in  this  situation.  For,  consider 

■•m  r_k 


g(t)=i  M[f(y)/f(t)]dy=  M[f(w  + t)/f(t)]dw=;  Mexpi  . 0.  [U  (w+  t)-U.  (t)]ldw. 
“C  ‘o  "0  l 1=1  1 1 1 J 

Therefore  g'  (t)  =J  M expj£^=1  0i[Ui(w  + t)-UL(t)  ]]^^=1  ©i  fU^(w+t ) -U[(t)  ]dw 
-fMejcpjzJ_10i[U1(y)  - U.(t)  ]}  [T,(t)  -H(y)]dy, 

and  (2.4)  holds,  with  st  (y)  = exp^E^^  9^ [U^(y)  -Ui(t)]|>  0.  We  have 
established  the  following  result. 

COROLLARY  1.  Suppose  f has  the  exponential  family  form  (2.7).  Then 
the  assertions  (a)  through  (d)  of  the  Theorem  hold,  where  T|'  is  given 
by  (2.8). 
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It  should  be  noted  that  the  form  (2.7)  is  quite  general.  Typically 
the  vector  0 = (©^,...,0^  is  a vector  of  parameters.  In  using  Corollary  1, 

. k 

a density  of  the  ostensibly  more  general  form  f (t)=C (0)h (t )exp[Z._ ^0. U . (t ) } 
should  be  written  as  f (t)=C (0*)exp(E  ^0*  U*  (t)},  where  k*  = k+1, 

8*  = » • • • »8k> 1) » and  U*(t)  = U^t),  i=l,...,k,  with  U*+1(t)=log  h(t). 

The  Theorem  is  applicable  also  to  certain  exponential  family  mixtures. 
Suppose  f can  be  expressed  as 


f (t)  = p f1(t)  + q f2(t), 


(2.9) 


where  0 < p < 1,  q = 1-p,  and 


fj(t)  = Cj(0)exp{E^=10ijU.(t)},  0 < t < M, 


(2.10) 


with  each  U^(t)  being  twice  differentiable  on  (0,M).  We  note  that 


Tl(t)  = f'(t)/f(t)  = - [f’(t)  + cf^(t)]/[f1(t)+cf2(t)],  where  c is 


defined  by  c = q/p.  Since  f^t)  + c f2(t)  = c C2(0)exp{E^10.2Ui(t)}?(t)  , 
where  §(t)  is  defined  by 


5(0  = 1+  [C1(0)/cC2(0)]exp{Ej=1(0.1  -0.2)U.(t)}, 


(2.11) 


it  follows  that  T|(t ) = - [?'(t)/§(t)  + (t)] . For  computational 

purposes  it  is  useful  to  note  that 


n(t)  = - (2j.1eilu[(t)  + [2^=1(0i2-0tl)u^(t)]/5(t)}, 


(2.12) 


We  establish  the  applicability  of  the  Theorem  by  showing  that  the  condition 


I’M 


(2.4)  always  holds  here.  Since  g(t)=  [f (w+t)/f (t) ]dw,  we  have 

-M  -0 


g’(t)=.  [f (w+t)/f (t)]'dw. 
°0 


But  f(w+t)/f(t) 


[f1(w+-t)+c  f2(w+t)]/[f1(t)+c  f2(t)]-[?  (w+t)/?(t)]exp{Z^=1012[Ui(w+t)-U1(t)]). 


Therefore  [f (wft)/f (t)] ,=exp{Z^=1012[U1(w+t)-U1(t) ]}{?' (wft)/?(t) 


§ (vH-t)|  '(t)/[§(t)]  +[§  (v+t)/§(t)  ]^i=i©t2  Consequently, 
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g'(t)='Mexp{^=1ei2[ui(y)-u.(t)]}[§(y)/§(t)]{?'(y)/?(y)-§,(t)/§(t) 

+ £ (y)-U' (t)]]dy  = s (y)[TKt)-TKy)]dy,  where 

1=  1 1 i ° t ^ 

st(y)=exp{E^=iei2[Ui(y)-Ui(t)]} [?(y)/?(t)]  > 0 for  all  t and  y in 
(0,M).  We  therefore  have  the  following  result. 

COROLLARY  2.  Suppose  f has  the  mixture  of  exponential  families  form 
described  by  (2.9)  and  (2.10).  Then  the  assertions  (a)  through  (d)  of 
the  Theorem  hold,  were  7]'  can  be  obtained  from  (2.12)  and  (2.11). 

In  concluding  this  section,  we  present  an  approach  which  can  speed 
implementation  of  the  Theorem  and  Corollaries  in  situations  (c)  or  (d), 
where  (2.5)  or  (2.6)  hold,  respectively.  Suppose  (2.5)  holds.  The 
crucial  issue  is  whether  (i)  or  not  (ii)  a number  (0,M)  exists 

satisfying  g'(yg)  = 0.  From  the  proof  of  (c),  it  is  clear  that  either 
g(t)  has  a single  extremum  in  (0,M),  namely  a maximum  (case  (i),  (B)), 
or  g(t)  is  strictly  monotone  decreasing  (case  (ii),  (I)).  The  crucial 
issue  can  therefore  be  stated  in  the  typically  more  tractable  form  of 
whether  (i)  g is  increasing  initially  (i.e.,  in  a neighborhood  of  0)  or 
(ii)  g is  decreasing  initially.  To  appreciate  the  usefulness  of  this 
representation,  consider  the  behavior  of  g near  0.  Suppose  limt^gg(t) 
exists  and  equals  0.  (This  is  equivalent  to  limt^gf(t)  = ®.)  Then  since 
g(t)  > 0 for  all  t€  (0,M),  we  conclude  that  g must  be  increasing 
initially,  (i)  holds,  and  f is  BT.  Similarly,  lii^gf(t)  = 0,  or 
limtjgg(t)  = ®,  implies  (ii)  holds,  and  f is  IFR.  If  f(t)  does  not 
tend  to  0 or  ®,  it  is  then  expedient  to  consider  the  behavior  of 
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g(t)Tl(t)  near  0.  Suppose  5 = limt^0g(t)T|(t)  exists,  possibly  equal  to 
-poo  or  . From  (2.2),  we  note  that  6 > 1 implies  g'  (t)  >0  in  a 
neighborhood  of  0,  so  that  (i)  holds,  and  f is  BT.  On  the  other 
hand,  5 < 1 implies  (ii)  holds,  and  f is  IFR.  If  6=1,  higher  de- 
rivatives of  g near  0 may  be  investigated  to  determine  whether  g 
is  initially  increasing  or  decreasing.  The  situation  of  (2.6)  is 
analogous  to  that  of  (2.5),  with  the  issue  being  whether  g is  de- 
creasing initially  (case  (i),  (U))  or  increasing  intially  (case  (ii) , 
(D)).  The  following  summarizes  our  modified  approach. 

LEMMA.  Suppose  (2.5)  or  (2.6)  holds  in  the  Theorem  or  in  Corollary  1 
or  2. 

(a)  Suppose  e = limt^f(t)  exists,  possibly  equal  to  0 or  ®>. 

(i)  If  e = 00  and  (2.5)  holds,  then  (B). 

(ii)  If  € = 0 and  (2.5)  holds,  then  (I). 

(iii)  If  e = 0 and  (2.6)  holds,  then  (U). 

(iv)  If  e = =»  and  (2.6)  holds,  then  (D). 

(b)  Suppose  6 = limtl0g(t)T](t)  exists,  possibly  equal  to  +»  or 

(i)  If  6 > 1 and  (2.5)  holds,  then  (B). 

(ii)  If  6 < 1 and  (2.5)  holds,  then  (I). 

(iii)  If  6 < 1 and  (2.6)  holds,  then  (U). 

(iv)  If  6 > 1 and  (2.6)  holds,  then  (D). 
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3 . EXAMPLES 

In  this  section  examples  are  given  which  illustrate  the  applica- 
bility of  Corollaries  1 and  2.  In  the  process  new  BT  models  are 
introduced. 


EXAMPLE  1.  A generalization  of  the  gamma  and  truncated  normal  densities. 
Consider  the  exponential  family  of  densities  of  the  form 

2 

f (t)  = C(a,p,y)exp{-  art  - pt  + y log  t}  , 0 < t < ®, 


where  the  natural  parameter  space  is  the  union  of  £ (ar,p , y) : - ® < a <®  , 
p > 0,  y > -1}  and  £(or,p,y):ar  > 0,  p=0,  y > -l}.  The  special  case  p=0 
gives  the  class  of  gamma  densities,  which  includes  the  subclass  of  ex- 
ponential densities  (y=0).  The  special  case  p > 0,  y=0  gives  the  class 
of  truncated  normal  densities.  To  apply  Corollary  1,  note  that  k=3,  M=«, 
(0^02, 03)  = (or,  p , y) , and  (U^t)  ,U2(t)  ,U3(t))  = (-t,  -t2,  log  t).  There- 
fore Tj ( t ) = - I^_^0JJ^(t)  = at  + 2p-y/t,  and  T)'  (t)  = 2p+y/t2. 

Case  1.  p > 0.  (a)  if  y > 0,  then  T|'  (t)  > 0 for  all  t > 0,  so  that  by 

Corollary  1(a),  (I)  holds.  (b)  On  the  other  hand,  if  y < 0,  then  (2.5) 

1/2 

holds  with  tg  = (-y/2p)  . Since  limtjgf(t)  = “,  by  Lemma  (a,i)  we 

have  that  (B)  holds. 


2 

Case  2.  p = 0.  (Gamma  density.)  Here  T|*  (t)  = y/t  . (a)  If  y > 0, 
then  by  Corollary  1(a),  (I)  holds,  (b)  If  y< 0,  then  by  Corollary  1(b), 
(D)  holds.  (c)  If  y = 0 (exponential  density),  then  r(t)  = f (t)/R(t) = a, 
a constant.  These  of  course  are  well-known  facts.  See  Barlow  and 


Proschan  (1975). 


^ 
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EXAMPLE  2.  A generalization  of  the  beta  density. 

Consider  the  exponential  family  of  densities  of  the  form 

f(t)  = C(a,p)exp{o  log (t/M)  + p log(l-t/M)},  0 < t < M, 

where  the  natural  parameter  space  is  {(o;,p):Q(  > -1,  p > — 1} . In  fact, 

C(a,p)  = M V(a+  1 )T(p  + 1)/T(a+  p + 2) . The  case  M=1  is  the  usual  beta 
density.  The  case  a = p = 0 is  the  uniform  density  over  (0,M). 

Corollary  1 is  applicable  here,  with  k=2,  = (a,p),  and 

(U^t),  (t ) ) = (log (t/M) , log (1-t/M) ) . Consequently, 

T)( t ) = p/(M-t)  - or/t,  and  TJ'  (t)  = [p  t2  + o(M-t)2]/t2(M-t)2. 

2 2 

Case  1.  a<  0,  p>0.  Define  the  quadratic  h(t)  by  h(t)  = pt  + ar(M-t) 

2 2 2 2 
= (or  + p)  t - 2M  a t + or  M . Note  that  h(0)  = a M < 0 and  h (M)  = p M >0. 

(a)  If  or+ p = 0,  then  h(t)  = -oM(2t-M),  so  that  (2.5)  holds  with 

tQ  = M/2.  (b)  If  or+p>0,  then  since  h has  its  minimum  at 

t = M a/ (a  + B ) < 0,  h is  strictly  increasing  in  (0,M),  and  crosses  0 

somewhere  in  (0,M).  Therefore  (2.5)  holds.  (c)  If  a+p<0,  then 

since  h has  its  maximum  at  t = Ma/(a+p)  > M,  h is  strictly 

increasing  in  (0,M),  and  crosses  0 somewhere  in  (0,M).  Therefore 

(2.5)  holds.  We  conclude  that  (2.5)  holds  for  any  (a,p)  satisfying 

a<  0,  p > 0.  Since  also  limt^f(t)  = ®,  we  deduce  from  Lemma  (a,i) 

that  (B)  holds. 

Case  2.  a>0,  p<0.  An  argument  analogous  to  that  of  Case  1 shows  that 

(U)  holds  here. 

Case  3.  a>0,  p > 0 or  a>0,  p>0.  From  Corollary  1(a)  it  is 

apparent  that  (I)  holds. 



•-  ■ . — J 
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Case  4.  a<0,  p<0  or  Qf<  0,  p<0.  From  Corollary  1(b)  it  is  apparent 

that  (D)  holds. 

Case  5.  a=  0,  p=0.  Here  r(t)  = f(t)/R(t)  = l/(M-t),  so  that  (I) 
holds . 

EXAMPLE  3.  The  lognormal  density. 

The  lognormal  density, 

_____  i I n 

f (t)  = (/2rr  a t)  exp{  - 2 [ (log  t-p.)/c]  },  0<t<®,  - ®<^<®,  o>0, 

2 

may  be  written  as  f(t)  = C(or,p)exp(-  a(log  t)  + plogt],  where 
2 2 

ar=  l/2a  >0  and  p = \i/ a -1,  -®<p<®.  Corollary  1 applies  with 

M=®,  k=2,  (ei,02>  = (°f,P)»  and  (U1(t)  ,U2(t))  = (-(log  t)2,  logt).  Thus 

2 

T]'  (t)  = (2or+  p - 2a  log  t)/t  , and  (2.6)  holds  with  t^  = exp{  1 + p/2or}  . Since 
limt , gf (t) = 0,  by  Lemma  (a,iii)  we  have  that  (U)  holds.  Although  this 
result  is  well-known  (see  Mann,  Schafer,  and  Singpurwalla  (1974)),  our 
approach  is  notably  swift. 

EXAMPLE  The  inverse  Gaussian  density. 

The  inverse  Gaussian  density, 

3 1/2  2 2 i 

f(t)  = (X/2rrt  ) exp(- X(t-p,)  /2p.  t],  0<t<®,  p,>0,  X>0, 

has  the  exponential  family  form  f (t ) = C (p,,X)e.xp{- (X/2p,  )t-(X/2)t  -(3/2)log  t}, 

which  is  amenable  to  Corollary  1 for  the  case  M=®,  k=3,  (e^.e^Q^)  = 
(X/2|i2,X/2,  - 3/2)  , and  (l^ (t) ,U2(t) .U^ (t ) ) = (- 1 , -t  1,logt).  Therefore 
T)*  (t)  = (X  - 3t/2)/t"^,  so  that  (2.6)  holds  with  tg=2X/3.  Since  lim^^f  (t)  = 0, 
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by  Lemma  (a,iii),  we  have  that  (U)  holds,  a fact  demonstrated  with 
greater  computational  effort  by  Chhlkara  and  Folks  (1977). 

EXAMPLE  5.  A cubic  exponential  family. 

Perhaps  the  most  straightforward  way  to  have  (2.5)  or  (2.6)  hold 
is  to  make  T) ’ linear.  This  is  accomplished  in  the  exponential  family 
case  by  the  class  of  densities, 

f (t)  = C(Qf,p,v)exp(-Qft  - pt  -yt  },  t > 0, 

where  the  natural  parameter  space  is  the  union  of  {(ar,p,y):  -®<a<®, 

-®<p<®,  y>0],  ((ck.P.y):  - ®<  a<  ®,  p>0,  y=0},  and  ( (a,p,y):a>  0, 

P = 0,  y=0}.  In  applying  Corollary  1,  we  note  that  M=®,  k=3, 

(ei,02,e3^  = and  (Uj(t:)»U2(t),U.j(t))  = Therefore 

2 

T)(t ) = ot+  2pt  + 3yt  , and  T]1  (t)  = 2p  + 6yt. 

Case  1.  p > 0,  y>0  or  P>0,  y>0.  Here  by  Corollary  1(a),  (I)  holds. 

Case  2.  p=Y=0-  As  noted  in  Example  1,  for  this,  the  exponential  density, 
we  have  constant  failure  rate  equal  to  a. 

Case  3.  p<0.  Then  (2.5)  holds  with  tQ=-p/3y.  Since  lim  gf(t)  = 

C(3«P>Y)  4 0,“,  we  cannot  apply  Lemma  (a).  We  therefore  consider 
6 = limtl0g(t)T|(t)  = or/C(a,p,y).  (a)  Suppose  a<  0.  Then  6<0<1,  and 
by  Lemma  (b,ii)  we  have  that  (I)  holds.  (b)  Suppose  a>0.  Then 
6=or  exp{-ort  - pt  - yt  }dt,  and  from  Lemma  (b,i)  and  (b,ii)  we  conclude 

v 0 

that  (B)  holds  if  6<1,  and  (I)  holds  if  6>1.  Further  insight  is  gained 

CO 

<*  2 3 2 

by  expressing  6 as  6 = exp(-y+  \y  -py  ]dy  = h(\,p),  where  X — p / or  >0 

« 0 


l 
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and  p=v/a  >0.  Therefore  (d/d\)h(\,p)  = y2exp{-y  + \y  - py  }>0,  so 

- 0 

that  h(\,p)  ts  strictly  increasing  with  \ for  fixed  p.  Since 

lim  h(X,p)=  exp{-y  - py  ]dy<  exp(-y)dy  = 1,  and  lim  h(\,p)  = ®, 

•>  o -0  A 00 

we  conclude  from  continuity  that,  holding  p fixed,  for  sufficiently  small 
\ we  have  6<1  (I),  and  for  sufficiently  large  \ we  have  6>1  (B). 
Similarly,  it  is  seen  that  h(A.,p)  is  strictly  decreasing  with  p for 
fixed  \,  and  that,  holding  \ fixed,  we  have  6<1  (I)  for  sufficiently 
large  p and  6>1  (B)  for  sufficiently  small  p.  The  case  6 = 1 is 
solvable.  Here  limt ^g 1 (t ) » 0.  However,  differentiation  of  (2.2)  gives 
limti0g"(t) = 2p/C(a,p  ,y)<0.  Thus  g is  initially  decreasing,  and 
(I)  holds. 


EXAMPLE  6.  Gamma  mixtures  with  common  scale  parameter. 

Consider  the  mixture  of  densities,  f (t ) = p f ^ (t ) + q f j(t ) , 0<p<l, 
q=  1-p,  where  f ^ , j = 1,2,  is  gamma  distributed  with  shape  parameter  Yj 


and  scale  parameter  or,  i.e.. 


r Yi 

f,(t)=  a J/r(v,)  t 


Yj  1 -at 


t > 0,  a>  0,  Yj>°* 


For  integer  valued  Yj « f^  is  the  density  of  the  time  until  the  occurence 
of  the  Yj-th  event  of  a Poisson  process  with  parameter  a.  If  an  item  fails 
upon  the  accumulation  of  precisely  Yj  Poisson  (a)  occurrences  (say  shocks), 
then  fj  is  the  appropriate  density  function  for  the  item's  lifetime.  If 
a population  consists  of  a mixture  having  proportion  p of  items  failing 
upon  Y^  Poisson  (a)  shocks  and  proportion  q»  1-p  of  items  failing  upon 


Yj  Poisson  (a)  shocks,  then  f is  the  appropriate  density.  We  generalize 
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here  to  all  positive,  not  necessarily  integer  valued  Yj  • We  further 
assume  without  loss  of  generality  that  y^>Y2« 

We  are  motivated  to  look  for  possible  BT  mixtures  by  considering 
the  following  representation  of  r(t).  Suppose  Y2<  1 and  y^>l,  i.e., 
f^  is  IFR  and  f^  is  DFR.  We  have 

r(t)-  [pf1(t)+q  f2(t)]/[PRl(t)  + q R2 (t)  ] - hj (t)^ (t)  + h2(t)r2 (t) , 

where  R^  is  the  reliability  and  r^  the  failure  rate  for  f , h^(t)=  [l+ch(t)l 
h2(t)  - l-h^(t),  and  h(t)  = R2(t)/R1(t),  with  c = q/p.  Thus  r(t)  is  a 
weighted  average  of  the  component  failure  rates,  where  the  weights  vary 
with  t.  In  fact,  it  is  easily  shown  that  h(t)  is  strictly  decreasing, 
ranging  from  1 at  t»0  to  0 at  t*».  Thus  the  weight  function  h^(t) 
increases  from  p to  1,  whereas  h2(t)  decreased  from  q to  0,  Since 
limtiQr2(t)  = «,  but  lim^^r^t)  = 0,  r(t)  has  an  initial  DFR  character. 

Since  lim  r.(t)«flr  and  lim  h.  (t) = 1,  r(t)  has  an  ultimate  IFR 
character.  We  shall  show  by  Corollary  2 that  indeed  f is  BT  for  this 
situation. 

Y. 

To  use  Corollary  2,  we  write  f as  f (t)»  a J/r(Y«)  exp{-at+(Y.-l)log  t] , 

J J «-  J J J 

t > 0,  which  gives  M=»®,  k-2,  (®11«®12»®21  ®22^  = (°«®*Y1-1.Y2‘1) » and 
(U^ (t) ,U2(t ) ) ■ (-t , log t)  . From  (2.11),  we  obtain  5(t) = 1 + b ta,  where 
a5Y^”Y2>0  and  b * (aa/c)r(Y2)/r(Y^)  > 0.  From  (2.12),  we  obtain 
T](t)  =•  a - e/t  + a/t  (1+b  t*)  , where  e*Y^"l*  Consequently, 

Tj*  (t)  - [t^(l+bta)^]  [eb^t^a+  [2e-a(a+l)]b  ta  + e-a].  In  investigating 

T|'  via  Corollary  2 it  suffices  to  consider  the  behavior  of  the  quadratic 

2 a 

h (w)  = ew  + [2e  - a(a+l)  ]w  + e - a,  where  w-w(t)  = bt  . 
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Case  1.  y^<1.  Here  e<0.  If  e=0,  then  h(w)  = -a(a+l)w-a<  0 for  all 

w>0.  If  e<0,  then  since  h(0)  = e-a<0  and  h'(w)<0  for  all  w>0, 
we  have  h(w)<0  for  all  w>0.  Therefore  T)' (t)<0  for  all  t>0  which 
implies  by  Corollary  2(b)  that  (D)  holds.  This  fact  is  consistent  with 
the  general  result  (see  Barlow  and  Proschan  (1975))  that  mixtures  of  DFR 
distributions  are  DFR. 

Case  2.  \^>1,  \^<  1.  Here  h(0)  = Y2_l<0»  which  implies  that  for  some 

Wq>0,  we  have  h(w)<0  for  all  w<Wq,  h(Wg)=0,  and  h(w)>0  for  all 
w>Wq.  Therefore  (2.5)  holds  with  tg  =*  (Wq/1j) Since  limtigf(t)  = «, 

it  follows  from  Lemma  (a,i)  that  (B)  holds.  To  summarize,  if  one  density 
is  strictly  IFR  and  the  other  is  strictly  DFR,  then  the  mixture  is  BT. 

Case  3.  Y^>1,  = 1.  Here  f^  is  the  exponential  density,  h(0)  = 0,  and 

h'(0)  = a(l-a).  If  a<  1,  then  h'(0)>0,  so  that  h(w)>0  for  all  w>0. 

If  a=l,  then  h'(0)=0.  However,  h"(w)  = 2e  implies  that  h(w)>0  for  all 
w>0.  Therefore,  if  a<l,  we  have  Tj'(t)>0  for  all  t>0,  so  that  by 
Corollary  2(a),  (I)  holds.  On  the  other  hand,  if  a>l,  then  h'(0)<0, 
which  implies  that  (2.5)  holds.  Since  lim^ ^f  (t ) * qar 4 0,  ®,  we  cannot 
use  Lemma  (a).  However,  since  6=  l/q>l,  I,emma  (b,i)  shows  that  (B)  holds. 
To  summarize,  if  f^  is  strictly  IFR  and  f^  is  exponential,  then  the  mixture 
is  IFR  if  V1  < 2 and  BT  if  Vx>2. 

Case  4.  Here  both  densities  are  strictly  IFR,  and  h(0)>0. 

Since  h'(w)=  2sv+2e  - a(afl)  is  an  increasing  function  of  w,  it  follows 
that  if  h'(0)>0,  then  h(w)>0  for  all  w>0.  Now  h'(0)>0  if  and 
only  if  d >a(a-l)/2,  where  d=e-a.  If  a<l,  this  inequality  is 
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satisfied  for  all  d.  Suppose  a>l  and  d<a(a-l)/2,  so  that 

h'(0)<0.  Therefore  h(w)>0  for  all  w>0  if  and  only  if  h(w)  = 0 

has  no  real  roots.  But  no  real  roots  exist  if  and  only  if 

2 2 

[2e-a(a+l)J  - 4e(e-a)<0,  i.e.,  (a-1)  <4d.  Therefore,  by  Corollary 

2(a),  we  have  that  (I)  holds  for  the  following  cases:  a<l;  a>l, 

, 2 

d>a(a-l)/2;  and  a>l,  (a-1)  /4<  d<  a(a-l)/2.  For  the  case  a>l, 

2 

d=  (a-1)  /4,  we  have  for  Wg=2a/(a+l)  that  h(w)>0  for  all  positive 

w^Wq,  with  h(Wg) = 0.  It  follows  that  for  tg  = (Wg/b)^a,  we  have 

TV  (t)>0  for  all  positive  ti<tg,  with  Tj*  (tg)  = 0.  From  (2.4)  it  then 

2 

follows  that  (I)  holds.  For  the  remaining  case  a>l,  d<(a-l)  /4,  it  is 
apparent  that  T|'  is  initially  positive,  then  negative,  and  ultimately 
positive.  The  results  of  Section  2 are  therefore  not  applicable.  However, 
we  may  rule  out  f as  BT  or  DFR,  since  lim^gf  (t)  = 0,  which  implies 
g(t)  is  initially  decreasing.  Further,  since  T)  is  ultimately  increasing, 
we  deduce  from  (2.4)  that  g is  ultimately  decreasing.  This  rules  out 
UBT.  A reasonable  conjecture  seems  to  be  that  f is  IFR  in  this  case 
as  well. 

EXAMPI£  7.  General  gamma  mixtures. 

We  now  generalize  the  mixture  of  Example  6 to  the  case  of  unequal 
scale  parameters,  i.e.,  for  j-1,2, 
r V,  -.  V.-l  -o  t 

fj(t)  = ^j  J/T(vj)  Jt  J e J , t>0,  Ofj  > 0,  Vj>0. 

The  analysis  for  this  situation  is  more  complicated,  and  we  restrict 
consideration  to  the  case  OCv^KY^®.  which  is  BT  in  the  « Or^ 
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context.  We  use  the  same  definitions  and  assignments  given  in  Example  6, 

r vi  V 

except  now  0^=0^,  and  b = ^al  ^ca2  y r(v2)/r(v1)  > 0.  Further, 

\ t fl 

we  define  X^cr^-0^*  Then  T](t)  = ar^-e/t  + (a+Xt)/t  (1+b  e t ) , so  that 

T|'(t)=  [t2(l+w)2]  (e(l+w)2-(a+Xt)2w-a(l+w) ] , where  w=w(t)=be*tt1, 

2 2 

Use  of  Corollary  2 involves  investigation  of  h(t) = e(l+w)  -(a+Xt)  w-a(l+w). 

Case  1.  X<0.  Here  w(0)  = w(<*>)  = 0,  and  h(0)  = h(“)  = £-a<  0.  Now  1) 

ultimately  decreasing  rules  out  IFR  and  BT,  and  limt  (t)  = 00  rules  out 
IFR  and  UBT.  DFR  remains  a possibility.  In  fact,  h(t)<0  for  all  t>0 
would  imply  (D)  from  Corollary  2(b).  However,  investigation  of  h at 
extrema  shows  that  h(t*)  = [l+w(t*) ] [e [l+w(t*) ] -a} , where  t*^-a/X.  Thus 
for  suitable  choices  of  (a,b,e),  h(t*)>0  and  Corollary  2 does  not  apply. 

Case  2.  X>0.  Here  w(0)  = 0 and  w(«)  =»  ®,  so  that  h(0)<0  and  h(®)>0. 

It  is  convenient  to  express  h(t)  as  h(t)=as(m),  where  m=Xt/a.  Then 

2 2 am  3 

s(m)  = d(l+z)  - a(m+l)  z-z-1,  where  d=e/a,  (0<d<l),  z=z(m)  = p,e  m , 

£ 

and  ^sb(a/X)  . Since  s(0)<0  and  s(«)>0,  there  exists  at  least  one 

positive  root  to  s(m)=0.  Suppose  is  such  a root.  In  order  to  show 

that  (2.5)  holds,  it  suffices  to  show  that  s' (nig)  > 0.  For  any  m, 

2 * 

s' (m)  = a[(»fl)/m]  [2d(l+z)z-2mz  - a(nH-l)  z-z].  Since  s(mQ)  = 0,  we  have 
2 2 

2dZg  = 2a(mQ+l)  zQ+ 2zQ+ 2 - 2d-4dzQ  (where  ZgSzOng)),  so  that 

s 1 (®q)  - a[(m0+l)/m0](z0[a(m0+l)2-2d  - 2mQ+l]  + 2(l-d)}.  Consider  the 

2 

quadratic  v(m)=a(nH-l)  -2d-2nH-l.  (a)  Suppose  a>l.  Then 

v(0)  =*  a-2d+l > 2(l-d) > 0,  and  v(m)  is  strictly  increasing  for  positive  m. 
Therefore  v(mg)>0,  which  implies  s'(mg)>0,  and  (2.5)  holds.  From 
Lenma  (a,i)  we  conclude  that  (B)  holds.  (b)  Now  suppose  a<l.  In 
this  situation,  unlike  the  a>  1 case,  the  number  of  positive  roots  to 


t 
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s (m)  = 0 may  exceed  1.  This  is  because  v' (0)  = 2(a-l) < 0 allows  for  a 
root  mQ  satisfying  s'(m0)<0.  As  an  example,  consider  a =3/16,  e=3/32, 

I*  = 2 exp(-3/16) , Then  s(l)=0  but  s'(l)  = -9/16.  In  fact  in  this 
example,  s(m)=0  has  3 positive  roots  so  that  T|'  is  initially  negative, 
then  positive,  then  negative,  and  finally  positive.  We  may  summarize  the 
case  a< 1 by  stating  that  (B)  holds  if  there  is  precisely  one  positive 
root  to  s(m)=0.  If  the  number  of  such  roots  exceeds  1,  our  general 
results  are  not  applicable,  although  it  is  easy  to  see  that  IFR,  DFR,  and 
UBT  are  ruled  out.  A seemingly  reasonable  conjecture  is  that  f is  BT 
here  also. 

We  conclude  this  example  by  considering  the  special  case  Y2 = 1 

(exponential  density),  with  a>l  and  X>0.  Recall  from  Example  6 

that  for  X = 0 , f is  BT  if  a>l  and  IFR  if  a=  1.  We  shall  discover 

that  in  the  present  case,  f is  BT  if  a>l  and  either  IFR  or  BT  if 

a=l.  Note  that  for  V2=*»  we  have  d=l.  Therefore  s(0)  = 0 and 

s(co)  = oo.  The  existence  of  a positive  root  to  s (m)  = 0 is  not  assured. 

If  one  exists,  say  m^m^,  the  argument  given  above  for  a>l  still 

applies  to  show  that  v(ii1q)  > 0 and  v'(mg)>0,  and  thus  that  (2.5)  holds. 

On  the  other  hand,  if  no  positive  root  exists,  we  conclude  that  h(t)>0 

for  all  t>  0,  so  that  (I)  holds.  To  show  existence  of  a positive  root,  it 

suffices  to  show  s is  initially  decreasing.  Since  z'  = z' (m) = a[(nH-l)/m]z, 

2 

we  may  express  s 1 (m)  as  s'(m)  = z'i|r,  where  i|i  = ijf(m)  s l+2z  - 2m  - a(nrt-l) 

Therefore  s' (0)  = a(l-a)  lim  ..(z/m)  =0,  since  lim  iri(z/m)=p,  if  a*  I, 

mtu  mlU 

= 0 if  a>l.  Consider  therefore  s"(0).  We  have  s"(m)  = z'  ijf 1 + z"i(f. 


and 
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with  z"  = z "(m)  = a [a+2a/m+  (a-l)/m  ]z.  (a)  Suppose  a=l.  Then  z' (0)  = p,, 

f (0)  = 2p,-4,  z"(0)  = 2p,,  and  |(0)  = 0,  so  that  s" (0) = 2p(p-2) < 0 only 

if  p,<2.  It  can  be  shown  that  s"'(0)=  12  for  p,=  2.  We  conclude  that 

for  a=  1,  s is  initially  decreasing  only  if  p,<2.  (b)  Suppose 

a > 1 . Then  z'(0)  = 0,  i|r' (0)  = -2(a+l)  , z"(0)  = a(a-l)  lim  _(z/m2),  and 

miu 

2 

i(r(0)  = 1-a.  Since  limmi0(z/m  ) = » if  l<a<2,  = p if  a=2,  and 

= 0 if  a > 2,  we  have  that  s"(0)=-«<0  if  l<a<2,  s"(0)  = -2p, < ( 

if  a=2,  and  s"(0)=0  if  a>2.  By  taking  higher  derivatives  it  can 

be  shown  that  for  a>n-l,  s^(0)  = 0 for  i<n-l,  and  s ^ (0)  = (l-a)z  (0)  , 

with  (0)  = a(a-l)  • • • (a-n+l)  lim  ^(z/m0) . Consequently,  for  all  a>l, 

we  have  s initially  decreasing.  In  all  instances  where  s initially 

decreases  (implying  (2.5)  holds),  we  need  to  consider  the  Lemma.  Since 

Y2  = 1>  we  have  lim^^f  (t)  = q a^,  and  Lemma  (a)  does  not  apply.  However, 

for  a>l,  limtjgT|(t)  = ar^  + \=  a^,  so  that  6=l/q>l,  and  Lemma  (b,i) 

concludes  (B).  For  a=l  with  p,<2,  lim^gTKt)  = + X-b  = c^b. 

2 

Therefore  6=  1/q  - p^^/qa^)  , which  may  be  less  than  or  greater  than  1, 

depending  on  (p,a^,a^) . The  following  summarizes  our  results  in  the 

exponential  case  (Y2=l)>  where  \>0  and  a>l.  If  a>l,  or  if  a = 1 , 

2 

p,<  2,  and  1/q  - p(a^/q  a^)  >1,  then  (B)  holds.  If  a=l  and  p,>2, 

2 

or  if  a=l,  p,<  2,  and  1/q  - p(a^/q  <1  , then  (I)  holds.  Since  the 
case  \=0  implies  p = ®,  these  results  are  consistent  with  the  \ « 0 
counterparts  in  Example  6. 

EXAMPLE  8.  Weibull  mixtures. 

In  Example  6,  we  showed  that  the  mixture  is  BT  for  a strictly  IFR 
gamma  density  and  a strictly  DFR  gamma  density  with  the  same  scale  parameter. 
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Here  we  show  that  the  corresponding  mixture  based  on  Weibull  densities 
is  not  BT.  This  disputes  a contention  of  Kao  (1959). 

Consider  the  mixture  f(t)=  p f^(t) +q  f (t),  where 


P ” 1 p 

fj(t)  = ap^(at)  j exp{-(at)  -*}, 


t > 0,  a>  0,  8 . > 0. 

J 


The  failure  rates  have  the  concise  form  r^  (t)  = a (at ) J . Now 
r (t)  = h1(t)r1(t)  + h2(t)r2(t)  , where  h1(t)=  [l+ch(t)]  1 and 

P P 

h2(t) = l-h1(t),  with  c = q/p  and  h(t) = R2(t)/R1(t) = exp{ (at)  1-(at)  2}. 

By  assuming  0<  P2<  1<  we  have  IFR  and  f2  DFR.  But  then 

limt,Qr(t)  = “ and  limt  ^r  (t)  = 0.  Therefore  f cannot  be  BT.  In  fact, 
it  is  seen  that  r(t)  initially  decreases,  then  increases , and  then  ulti- 
mately decreases  for  the  given  mixture.  We  may  generalize  to  unequal 
scale  parameters,  a^.  Again  lim^^rtt)  =®  and  limt  #r  (t)  = 0, 
which  rules  out  BT.  The  same  conclusion  holds  when  f2  is  Weibull  with 
a positive  valued  location  parameter.  Therefore  Kao's  contention  that 
such  Weibull  mixtures  are  BT  is  incorrect. 


4.  STATISTICAL  INFERENCE 

The  models  in  Examples  1 through  5 have  the  exponential  family  of 
densities  form.  As  a result,  statistical  inference  is  simplified  by  a 
reduction  to  sufficient  statistics,  and  optimum  hypothesis  tests  and 
confidence  region  procedures  can  be  obtained.  See  Lehmann(1959) . We 
consider  in  this  section  certain  point  estimation  techniques,  particularly 
for  the  model  of  Example  1. 
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THE  GENERALIZED  GAMMA -TRUNCATED  NORMAL  DENSITY  OF  EXAMPLE  1. 

The  parametrization  given  in  Section  3 obscures  the  role  of  ot  in  the 
(3  > 0 situation.  Since  |a|  is  a scale  parameter  unless  a=  0,  it  seems 
natural  to  break  the  overall  model  into  the  following  four  classes, 
which  effectively  delineate  these  roles  of  a.  l:(a>0,  0>O,  y>-l). 

2:  (a<  0,  p>0,  y>-l).  3:(d=0,  p>0,  y>-l).  4;(a>0,  p=0,  y>-l). 

For  convenience  in  estimation,  we  re-parametrize  slightly. 

Class  1:  f(t)  = [a^/r(9,p)]exp{-at  - 9(at)^}tp  \ or>0,  9>0,  p>0. 

f* 00  2 1 

F(9,p)  is  defined  by  r(9,p)=  exp[-y  - 9y  }yH  dy. 

" 0 

Class  2:  f (t)  = [ctP/ A (9,  p)  ] exp{at  - 9(ct  )2}t“  ^ , c>0,  9>0,  p>0. 

1 ® 2 Q — 1 

A(9,p)  is  defined  by  A(9,p)=  ^exp[y  - 9y  }yp  dy. 

Class  3:  f(t)=  [X?/H(p) ]exp{- (\t )2} tp  *,  X>0,  p>0. 

H(p)  is  defined  by  H(p)=  exp{-y  }yP  dy. 

* 0 

Class  4:  f (t)  = [orP/T(p)]e  attp  ar>  0,  p>0. 

Note  that  in  Classes  1,  2,  and  3,  f is  BT  if  p<l,  IFR  if 
p>  1,  and  truncated  normal  if  p=  1.  Class  4 is  the  class  of  gamma 
densities,  where  f is  DFR  if  p<I,  IFR  if  P>1,  and  exponential 
if  p=l. 

Maximum  likelihood  estimation  of  the  parameters  is  tractable  for 
this  model.  In  practice,  the  statistician  would  form  a collection  of 
possible  lifetime  models  by  choosing  some  or  all  of  the  four  classes. 
Typical  choices  might  be  Classes  [1,2,3, 4},  or  [l,2],  or  [l,4],  or  [2,4], 
or  {1,2} . For  each  class  included,  the  MLE  of  the  corresponding  parameters 
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would  be  computed,  as  well  as  the  value  of  the  likelihood  at  this  point. 

The  MLE  of  (a,p,Y)  would  then  be  the  value  of  (a,(3,Y)  corresponding  to 

the  MLE  of  the  parameters  in  the  class  with  the  largest  maximum  likelihood. 

For  example,  suppose  we  wish  to  estimate  (a,{5,Y),  assuming  the  collection 

of  possible  lifetime  densities  is  the  set  of  Classes  (l,4}.  Maximum 

likelihood  estimation  for  the  individual  classes  yields  the  MLEs,  say 

(a^,0^,p^)  for  Class  1 and  (o^, p^)  for  Class  4.  If  the  likelihood 

evaluated  for  Class  1 at  (a^.0,,p^)  exceeds  the  likelihood  evaluated  for 

2 

Class  4 at  (ar^,p^),  then  the  MLE  of  (a,{3,Y)  is  (0^,9^  Q^.p^-l) . If  the 

likelihood  (4)  at  (a, ,P.)  exceeds  the  likelihood  (1)  at  (a, ,9. ,p, ) , 

4 4 111 

then  the  MLE  of  (a,B , y)  is  (a, ,0,p  -1). 

*+  4 

We  present  now  the  individual  MLEs  for  the  four  classes.  Assume  a 

random  sample  of  failure  times  is  taken,  with  observed  values  t,,...,t  . 

1 n 

Computation  of  MLEs  is  simplified  by  the  exponential  family  form.  In 
general,  if 

f(t)  = C(0)h(t)exp{l^=10.U.(t)},  0 < t < M , (4.1) 

then  it  is  easily  seen  that  any  solution  to  the  set  of  equations, 

(3/39^) log  L = 0 , i=l,...,k,  (where  the  likelihood  L = II?_^f (t^ ) ) , 
is  a solution  to  the  set  of  equations, 

E Ut(T)  = St(t),  i=l,...,k,  (4.2) 

where  T denotes  the  time-to-failure  random  variable,  and 

( S^(t) S^(t))  is  the  sufficient  statistic  given  by  S^(t)  = ^ E?_jU^(tj). 
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Class  1.  Here  ( S1(t) ,S2(t) ,S3(t))  = (£  St j Zt2,^  2 log  t^)  , and 
( E U1(T),E  U2(T),E  U3(T))  = ( E T,  E T2,  E log  T)  . Now 

et  = [r(9,p+i)/or(e,p)l [ap+1/r(e,p+i) ]exp{-at-e(at)2}tpdt=r(e,p+i)/ar(e,p)„ 

"0 

2 2 

Similarly,  ET  = r(Q,p+2)/a  F(0,p)  . In  a like  manner,  E log  T = \jf  (e , p)-log  a , 
where  i|r(0,p)  s fT(0,  p)  ] ; exp{-y-0y  }yP  log  y dy  = (d/dp)log  T (0,p). 

c U 

We  may  solve  the  system  (4,2)  to  obtain  the  MLE  (a^,9^,p^).  A computational 

2 2 

version  of  the  system  is  the  set  of  equations  S^/S0=r  (0^ , p^+1  )/T(0^ , p^)T(0^ , p^+2] 

s3  - log  s1=  , pL>  - iog[r(01,p1+i)/r(01,p1)],  and  = r(01,p1+i)/s1r<01,p1). 

The  first  two  equations  may  be  solved  simultaneously  (by  computer)  to 
yield  (0^,p^).  The  third  equation  then  provides  a^. 

Class  2.  Here  the  S.(t)  and  U.(T)  are  identical  to  those  in  Class  1.  It 

i ~ i 

0 o 

is  easily  seen  that  E T = A (0,  p+l)/cA(0 , p)  , E T = A (0,  p+2)/a  A(0,  p) , and 

E log  T = £(0,p)-log  O,  where  £(0,  p)  5 [ A (0 , p ) ] exp{y-0y  }yp”  log  y dy 

u0 

= (S/Sp)logA(0,p) . The  resultant  computational  version  of  the  system 
(4.2)  which  yields  the  MLE  (c2>02, is  t2le  set  e<luations 

S1/S2=A2  (02*fi2+1)/A(92,P2)A(92’p2+2),S3‘logSl=C(e2*P2)‘log[A(e2,P2+IJ/A(e2’P2) 
and  a2  = A (02 , P2+l) /SL A (02 , P2) . 

Class  3.  Here  ( S.  (t),S  (t ))=(-;  I t2  ^ Z log  t ) , and  ( E U.  (T),E  U (T)) 
i ~ z ~ n j n j i i, 

= (ET2,ElogT).  Then  ET2  = H(p+2)/\2  H(p),  and  E log  T = W(p)-log  where 

1 * 2 P “ 1 

W (p)3[H(p) ] exp{ -y  }yP  log y dy  = (d/Sp)logH(p)  . A computational  version 
J 0 

of  (4.2)  yielding  the  MI£  (\3>p3)  is  the  set  of  equations  S3  - log  S2 
= W(p3)--|  log(H(p3+2yH(P3)]  , and  \2 = H(p3+2)/S2H (p3) . 
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Class  4.  Here  ( SL(jt)  ,S2(jt))  = St  E log  t^)  , and  ( E Uj_ (T)  ,E  U2(T>) 

= ( E T,E  log  t)  . Since  ET=  p/a,  and  E log  T = ty(p)-log  or,  where 

•^r(p)  = (9/Sp)l°g  T(p) , we  obtain  the  well-known  result  that  the  MLE 

is  determined  by  S^-log  = |(p^)-log  p^,  and 

Large  sample  confidence  regions  for  the  parameters  of  each 

individual  class  can  be  obtained  from  asymptotic  normal  theory.  Let 

0 denote  the  MLE  of  0,  based  on  a random  sample  of  size  n taken 

from  (4.1).  Under  suitable  regularity  conditions  (see  Zacks  (1971)), 

the  limiting  distribution  of  /n(0  - 0)  is  k-dimensional  multivariate 

normal  with  mean  vector  0 and  covariance  matrix  I ^(9),  where 
2 

1(0)  = - E(3  /d0t$0j)log  f.  It  is  straightforward  to  show  that  for 

Class  1, 


l(or,0,p)  = 


JL  + r_(e,P+_2) 

22  r(0,p) 

a a v * ' 


2 r(8,p+2) 

a r (9, p) 


1 

a 


2 F(9,p+2) 

a r(0,p) 


1 

a 


r(e,p) 


n^©i£+2)' 
L r(9»P)_ 


2 

[t(e.P)-*(e,<H-2)j  iog  r(©,p) 


ap 


(4.3) 


For  Class  2,  l(c>0>p)  is  identical  to  (4.3)  if  we  substitute  o for 
and  the  function  A for  f.  For  Class  3, 


j 
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I(X,P) 


-£-4.-1  HQH-2) 


2 + 2 

X X 


H(p) 


— 2 log  H(p) 
op 


As  is  well-known,  for  Class  4, 


I(a.P) 


2 

a 


J. 

a 


i 

a 


.2 

-2-5  log  r<p) 

sp 


THE  OTHER  EXAMPLES.  We  consider  briefly  questions  of  inference  for 
Examples  2 through  7 from  Section  3. 

In  Example  2,  if  T denotes  the  time-to-failure  random  variable, 
then  T/M  has  the  usual  two  parameter  beta  density  over  (0,1).  Assuming 
M to  be  known,  we  may  employ  standard  inference  procedures  for  (a,j}). 
See  Johnson  and  Kotz  (1970).  Available  procedures  likewise  handle  the 
lognormal  (see  Mann,  Schafer,  and  Singpurwalla  (1974))  and  inverse 
Gaussian  densities  (see  Chhikara  and  Folks  (1977)). 

For  the  cubic  exponential  family  case  of  Example  5,  we  see  from 
(4.2)  that  the  MLE  of  (o,(3,y)  is  the  method  of  moments  estimator. 
Computation  would  require  the  use  of  a computer,  due  to  the  complicated 
nature  of  C(a,(5,y). 
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For  Examples  6 and  7,  the  lack  of  non-trivial  sufficient  statistics 
makes  the  method  of  maximum  likelihood  forbidding.  Method  of  moments 
estimators  are,  however,  computable.  The  i-th  moment,  p,^,  of  f is 
= PM'li  + *1^21’  where  p,^  = r(Yj  + i)/or*T(Yj  ) Is  the  i-th  moment  of 
fj,  j-1,2.  The  system  of  equations,  p^  = ^ l"=1  t*  , i=l,...,k,  yields 

the  desired  estimates.  In  Example  6,  we  use  k=4  to  estimate  (a,Y^,Y2>P)» 
and  in  Example  7,  we  use  k=5  for  ,n2 , Yj^ , Y2 »P)  • Solution,  of  course, 
requires  use  of  a computer. 
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